Incoherent Bi off-centering in Bi2Ti2060' and Bi2Ru2060': Insulator versus metal 
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In the cubic, stoichiometric oxide compounds Bi2Ti2060^ (also written Bi2Ti207) and 
Bi2Ru2060' (also written Bi2Ru207) Bi^+ ions on the pyrochlore A site display a propensity to 
off-center. Unlike Bi2Ti2060^ Bi2Ru2060^ is a metal, so it is of interest to ask whether conduction 
electrons and/or involvement of Bi 6s states at the Fermi energy influence Bi^^ displacements. The 
Bi^^ off-centering in Bi2Ti2060^ has previously been revealed to be incoherent from detailed by 
reverse Monte Carlo analysis of total neutron scattering. Similar analysis of Bi2Ru2060^ reveals 
incoherent off-centering as well, but of smaller magnitude and with distinctly different orientational 
preference. Analysis of the distributions of metal to oxygen distances presented suggests that Bi in 
both compounds is entirely Bi"^^. Disorder in Bi2Ti206 0^ has the effect of stabilizing valence while 
simultaneously satisfying the steric constraint imposed by the presence of the lone pair of electrons. 
In Bi2Ru206 0^ off-centering is not required to satisfy valence, and seems to be driven by the lone 
pair. Decreased volume of the lone pair may be a result of partial screening by conduction electrons. 



PACS numbers: 61.05.fm, 61.43.Bn, 71.30.+h 

INTRODUCTION 

The oxide pyrochlores ^2^2060^ usually abbreviated 
^42^207, are well known for their ability to accommo- 
date magnetic cations on interpenetrating sub-lattices of 
corner-connected 0^744 tetrahedra and ^Oe octahedra. 
The geometry of the sublattices often result in magneti- 
cally frustrated ground states, [1, 2 resulting in spin glass, 
spin ice, or spin liquid phases. [3 Frustration of concerted 
atomic displacements (dipolar frustration) has also been 
proposed, when electronic dipoles rather than magnetic 
spins are placed on the pyrochlore A site, [4 -6 with the 
suggested "charge-ice" [7 displaying the appropriate en- 
tropic signatures. [8j Large atomic displacement param- 
eters associated with the A-site cation, seen in Fourier 
maps of La2Zr207[9^ and in Bi2M207 {B Ti, Zn, Nb, 
Ru, Sn, Hf), p!QUT3] contribute to the expanding body of 
evidence that pyrochlores prefer to accommodate cation 
off-centering via incoherent disorder rather than in or- 
dered non-cubic ground states. 

Among 612^207 pyrochlores, B = Ti, (Zn/Nb), Sn, 
and Hf are insulators, while B = Ru, Rh, Ir, and Pt are 
metals. [14] The behavior of Bi^+ is suggested to be quite 
different in insulating and metallic pyrochlores. In in- 
sulating Bi2Ti2060' and Bi2Sn207 (cubic above 920 K) 
the Bi is offset by ~0.4 A from the ideal site. This Bi 
off-centering is aperiodic, but otherwise analogous to the 
correlated motion of lone-pair active Bi^+ in the ionic 



conductor Bi2 03 [T5] or multiferroic BiFeOs. [16^ In 
metallic pyrochlores, the suggestion is that weaker A- 
O' interactions preclude any displacement at alL pTHTO] 
However, Rietveld refinements for compounds where B = 
Ru, Rh, and Ir show that Bi off-centering is still present 
in experiment. p!0^, 20-22 The short-range correlation of 
Bi displacements has been probed using reverse Monte 
Carlo modeling of the diffuse streaks in electron diffrac- 
tion patterns. [23 

In this contribution, we compare structural details 
in insulating Bi2Ti2060' and metallic Bi2Ru2 060' us- 
ing pair distribution function (PDF) analysis with least- 
squares and reverse Monte Carlo (RMC) modeling. 
These techniques reveal the precise structural tenden- 
cies of Bi^+ off-centering, even in the case of incoher- 
ent ice-like disorder. [6 Specific details of bond distances, 
angles, and real-space shapes can be extracted from the 
RMC model because it is predicated on fits to the atom- 
atom distances in the PDF. The work should be placed in 
the context of ionic off-centering in metallic systems and 
the screening of ferroelectric dipoles, as first suggested 
by Anderson and Blount, [25 , which finds application in 
the context of heavily-doped perovskite titanates.[26] We 
find that the large displacements in Bi2Ti2060' can be 
reconciled with bond valence analysis, but displacement 
in Bi2Ru2060' are not driven by valence considerations 
alone. 

To introduce the comparison, Fig.[l] shows that the 
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FIG. 1. (Color online) Heat capacities of three pyrochlores 
show a strong dependence on the magnitude of static A-siie 
disorder (large in Bi2Ti206 0' present in Bi2Ru2 06 0^ none in 
¥2X1207). Note that the electronic contribution to the heat 
capacity of Bi2Ru2 06 0^ has been subtracted. The results of 
Rietveld refinement of Bragg neutron scattering displayed as 
95% ellipsoids on (orange) and Bi (black) on Bi2Ru206 0^ 
(left) and Bi2Ti2060^ (right) point to static Bi^^ disorder 
being enveloped within the large disks. Data adapted from 
refs.glElEl]. 



scaled heat capacities of Bi2Ti2060' and Bi2Ru2060' 
display pronounced low-temperature humps in plots of 
C/T^ vs. T that are indicative of local Einstein modes, 
suggestive of glassy disorder. This large local-mode con- 
tribution is largely absent in Y2Ti207, which has no 
lone pair and no experimentally resolvable displacive 
disorder. [8 Disorder on the 0'Bi4 sub-lattice is also seen 
in the average structure Rietveld refinement of Bragg 
neutron scattering, masquerading as very large atomic 
displacement parameters (ADPs) on the Bi sites, dis- 
played as 95% ellipsoids in Fig.[T] 



METHODS 

Preparation and average structure analysis of 
Bi2Ti2060' powder used in this study has been reported 
by Hector and Wiggin.[llj Bi2Ru2060' was prepared as 
single crystals by Tachibana [22^ and finely ground prior 
to measurement. Time-of- flight (TOP) neutron powder 
diffraction was collected at the NPDF instrument at 
Los Alamos National Laboratory at 298 K and 14 K. 
Rietveld refinement made use of the GSAS code. [27]. 
Extraction of the PDP with PDFGetN |28| used Qmax 
= 35 A~^. Reverse Monte Carlo simulations were run 
using RMCProfile version 6. [29] Electronic densities 
of states (DOS) were calculated by the linear muffin-tin 
orbital method within the atomic sphere approximation 
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FIG. 2. (Color online) The computed total densities of states 
(DOS) (a,b) for Bi2Ti2060' and Bi2Ru2060'. Below, Bi s 
(dark) and p (light) states are shown. Electron localization 
functions (ELFs) for both compounds are shown alongside as 
gold lobes around black Bi ions. Both ideal structures show 
annuli of lone pairs around Bi ions. 



using version 47C of the Stuttgart TB-LMTO-ASA 
program. [30 Bond valence sums (BVS) were extracted 
from the RMC supercell as described in previous work 
on CuMn204,[3r using the Rq values of Brese and 
0'Keeffe.[32^ 



RESULTS AND DISCUSSION 

The total DOS for Bi2Ti2060' and Bi2Ru2060' are 
shown in Pig.|2|a,b). The features are similar, with 
metallic Bi2Ru2060' shifted in a nearly rigid-band fash- 
ion by approximately 2 eV downward, in agreement with 
previous work. [18 Partial Bi s and p DOS are shown in 
Fig.[2|b,c). In both cases, some Bi s states are present at 
the top of the filled Bi p and B d bands. Bi s states are 
plotted as electron localization function (ELF ^ 0.65) iso- 
surfaces at the right of Fig.[TJa,b). With Bi in their ideal 
16c positions, the ELFs both show essentially identical 
circularly averaged lone pairs. [7 Assuming similar lone 
pair-cation distances, the ELFs imply that both com- 
pounds should have the same cation displacement in the 
real ground state structures. 

Time-of-flight neutron diffraction Rietveld refinements 
at T = 14 K and 300 K were performed using the ideal 
pyrochlore model with Bi on the 16c sites and anisotropic 
ADPs. Fits at 14 K are shown in Fig.|3] No substantial 
differences were found from the analysis of Hector and 
WigginpT] or Tachibana et a/. [22 including the occu- 
pancy; stoichiometric Bi2Ti2060' and a Bi occupancy 
of 0.97 for Bi2Ru2060'. The ADPs from 14 K Rietveld 
refinement are displayed as 95% ellipsoids in Fig.[T] The 
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FIG. 3. (Color online) Time-of- flight neutron diffraction 
Rietveld refinement of Bi2Ti2060^ (top) and Bi2Ru206 0^ 
(bottom) at 14 K using the ideal pyrochlore structure with 
anisotropic thermal parameters. The fit to Bi2Ti2 060^ is 
visibly worse due to large amounts of diffuse shoulder inten- 
sity. This diffuse scattering is a result of large local Bi and 
displacements. 



most apparent difference is the larger, disk-shaped ellip- 
soid representing Bi in Bi2Ti2 060^ Bi is known to be 
off-centered from the ideal 16c position in a ring normal 
to the O'-Bi-O' bond. [6] HTJ [33] The large anisotropic 
ADPs of Bi envelop this ring. A split Bi position can give 
a better fit to the diffraction data. The previous study 
found that there is a slight tendency for Bi to prefer the 
96h positions. [6 This represents a six- fold splitting of the 
Bi into sites that are displaced ~ 0.4 A from the ideal 
site, and pointing between nearby O ions on 48/ sites. Bi 
ADPs in Bi2Ru2060' also suggest displacive disorder. 
They too are anisotropic and appear as slightly flattened 
ellipses in Fig.[l] 

The most straightforward way to compare the 
Rietveld-refined unit cell with the local structure is via 
least-squares PDF refinements, shown in Fig.[4j Two is- 
sues should be considered. First, the fit for Bi2Ti2060' 
is significantly worse overall than Bi2Ru2060^ This 
implies that the local structure of Bi2Ti2060' is more 
poorly described by the Fd3m unit cell. Second, the fit 
for Bi2Ti2060' is poorest at low r, which contains details 
of near-neighbor atomic distances (of particular impor- 
tance are Bi-0 and Bi-0'), and is still unsatisfactory at 
higher r even though a larger number of pairs are being 
included. The fit for Bi2Ru2060' is equally decent at all 
r values up to 20 A. 

These PDF fits do not give the positions of atoms in 
either compound, but they quickly reveal valuable in- 
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FIG. 4. (Color online) Least-squares PDF fits (using the unit 
cell from Rietveld refinement) do not reproduce the low-r 
structure of Bi2Ti2 06 0^ due to the inability of a unit-cell 
based description to accommodate incoherent static displace- 
ments. The fit to Bi2Ru2060^ is significantly better because 
the static displacements are smaller. 

formation about the presence of local atomic displace- 
ments (more apparent in Bi2Ti2060' than Bi2Ru2 060') 
and the correlations between them (still unable to be 
averaged for r < 20 A). Extracting structural tenden- 
cies of geometrically frustrated compounds via least- 
squares refinement is inherently difficult because there 
is no straightforward way to model the large, complex 
collection of discrete displacements needed to reproduce 
the disorder. Least-squares is not a suitable algorithm for 
determining so many free Bi positions, especially when 
their interactions may be correlated. Instead, we remove 
symmetry constraints and use RMC to investigate how 
Bi are distributed within a large supercell. 

Simulations were carried out using the RMC method 
to investigate the precise positions of Bi. Simultane- 
ous fits to the PDF and Bragg profile for Bi2Ru2060' 
are shown in Fig.jsJ Unit-cell based modeling (least- 
squares refinements, including Rietveld) usually fails to 
model incoherent static displacements. In contrast, each 
RMC supercell contains thousands of ions of each type. 
Folding the RMC supercell into a single unit cell pro- 
duces "point clouds" of ions at each crystallographic site. 
These clouds display the propensity of ions to displace 
from their ideal positions. The mean squared displace- 
ment of points are in quantitative agreement with the 
average ADPs obtained from Rietveld refinement. Map- 
ping these points as two-dimensional histograms (Fig.[6| 
shows the tendency of Bi nuclei to offset in Bi2Ti2060' 
and Bi2Ru2060^ The Bi clouds are viewed normal to 
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FIG. 5. (Color online) RMC fits to the PDF and Bragg profile 
of Bi2Ru206 0^ at 14 K constrain the local structure (atom- 
atom distances) and long-range periodicity (cell parameter, 
ADPs) 



the O'-Bi-O' bond (top) and perpendicular (bottom) for 
two temperatures. 

The Bi ring in Bi2Ti2060' is evident from Fig. [g] and 
has a diameter of ~0.8 A. At 300 K, the ring is more 
diffuse. We attribute this to thermal broadening. Inter- 
estingly, the ring is not a perfect circle. It has a sixfold 
symmetry corresponding to the preference for Bi to oc- 
cupy the 96h positions (corners of the hexagon), which 
point between the six neighboring 48/ O ions in the TiOe 
network. [6] 

The Bi distribution in Bi2Ru2060' is distinctly differ- 
ent, but static displacement is still present. The displace- 
ments are densely clustered close to the ideal position and 
there is no hollow center as in Bi2Ti2060^ However, the 
perpendicular view reveals that the Bi distribution is still 
disk-shaped. Most surprising is the symmetry of the disk. 
It also has a hexagonal shape but the hexagon is rotated 
30° with respect to what is seen in Bi2Ti2060^ with flat 
edges on the left and right, and corners on the top and 
bottom. This implies that Bi is displacing toward the 
nearby 48/ O onto 96^ positions. The Bi offset roughly 
agrees with the value of 0.16 A found in the split-site 
model of Avdeev. [10 

Quantitative RMC Bi nuclear density as a function of 
the angle around a ring normal to the O'-Bi-O' bonds 
is shown in Fig.[6jb,c). Both compounds show sixfold 
modulation fit by a cosine curve with a period of 60°, 
but their oscillations are offset by 30°. 

Bond valence sums are calculated for each individ- 
ual cation in the supercell and plotted as histograms in 
Fig.[7|for (a) Bi2Ti2060' and (b) Bi2Ru2060' at T = 
14 and 300 K. In all cases, BVS distributions are cen- 
tered on the expected valence: Bi^+, Ti^+, Ru^+, and 
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FIG. 6. (Color online) Clouds of Bi nuclear intensity (a) 
for Bi2Ti2 06 0^ and Bi2Ru206 0^ are viewed along (top) and 
normal to (bottom) the 0^-Bi-O^ bond. Static disorder in 
both produces hexagonal ring or disk shapes centered on the 
ideal position. B-siie cations Ti and Ru are (bottom) are 
shown for size comparison. Distribution of Ti and Ru as a 
function of angle around the ring is shown in (b) and (c). 
The sixfold modulations indicate a preference for 96h and 96g 
sites. Data in (b) adapted from ref. [6]. 



0^~. A slight sharpening is seen for the low-temperature 
measurement. These distributions reveal that the RMC 
simulations contain chemically reasonable bond lengths 
despite the absence of such distance constraints in the 
simulations. They also reveal that there is no tendency 
for Bi^+ in Bi2Ru2060^ supporting the conclusion from 
LMTO calculations that Bi 5s states are localized far 
below the Fermi energy. 

The calculated Bi BVS versus displacement in Fig.[7|^c) 
demonstrates why displacements are more pronounced in 
Bi2Ti2060^ In the average structures, Bi2Ti2060' and 
Bi2Ru2060' respectively obtain only 1.31+ and 1.38+ 
per Bi from bonds to 0^ The majority of the valence is 
obtained from 48/ O and increases with Bi displacement, 
represented by the upward curve. Bi in both compounds 
gain about the same valence from 48/ O but a large dis- 
placement is required to reach the retracted Ti206 sub- 
lattice. The result is that each Bi in Bi2Ti2 060' gains an 
uneven amount of charge from the six 48/ O, locking in 
dipoles. This explains why the Bi displacements appear 
to be static in variable temperature measurements. [6] 
The rigidity of the ^206 sublattice is confirmed by small 
ADPs on 48/ O in both compounds, even in the presence 
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FIG. 7. (Color online) Bond valence sums (BVS) ex- 
tracted from the RMC super-cells show entirely Bi^^ in (a) 
BisTisOeO' and (b) BisRusOeO'. The small shoulder of 
is a consequence of peak broadening. Only Ti^^ (note 
peak center) exists in the sample. In (c), Bi BVS is plot- 
ted as a function of displacement into 96^ and 96h positions. 
Bi2Ti206 0^ requires significant off-centering to obtain valence 
from 48/ O anions. Dashed lines show the actual displace- 
ment and BVS. 



of nearby Bi off centering. In Bi2Ru2 060^ the RU2O6 
network pushes 48/ O closer to Bi so that valence is sat- 
isfied. 

We have considered whether covalency could lead to 
difficulties in using BVS to judge valence in Bi2Ru2060^ 
should more covalent bonding (shorter M-0 bonds) lead 
to Bi and Ru requiring more than the formal 3+ and 4+ 
to be satisfied? This seems unlikely, not only because 
Rietveld refinement and RMC find the desired states 
centered near the nominal values. The average struc- 
ture of semiconducting BiCaRu207 ^ displays much 
higher BVS sums (3.25+ for Bi, 4.20+ for Ru) than 
Bi2Ru2060^ but large anisotropic ADPs on the A site 
portend static disorder nonetheless. 

In Bi2Ti2060^ off-centering helps satisfy Bi valence 
and the lone pair can be accommodated in the oppo- 
site direction. In Bi2Ru2 060^ no off-centering is neces- 
sary from valence considerations, so static disorder may 
be driven by lone-pair activity. Metallic screening in 
Bi2Ru2060' is expected to decrease repulsions of the lone 
pair from nearby 0,[21 allowing Bi to stay closer to the 
ideal position than the traditional cation-lone pair dis- 
tance would dictate. This is supported by the idea that 
lone pairs often exhibit decreased volumes. |35j 

CONCLUSIONS 

In conclusion, reverse Monte Carlo stuctural analysis 
using total neutron scattering provides a detailed view of 



the incoherent static displacements in Bi2Ti2 060' and 
Bi2Ru2060^ Real-space maps of static displacements 
reveal the distinct magnitudes and directions of Bi off- 
centering in Bi2Ti2060' and Bi2Ru2060^ While static 
displacements in the insulator Bi2Ti2 060' can be under- 
stood on the basis of valence satisfaction alone, the cause 
for displacements in metallic Bi2Ru2060' is not captured 
by first-principles calculations on the ideal compound or 
by the bond valence sum. An incoherent lone-pair driven 
distortion is present but is partially screened by the con- 
duction electrons. 
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